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Abstract 



I review recent works on the problem of inducing large- N QCD by matrix fields. 
In the first part of the talk I describe the matrix models which induce large- N QCD 
and present the results of studies of their phase structure by the standard lattice 
technology (in particular, by the mean field method). The second part is devoted 
to the exact solution of these models in the strong coupling region by means of the 
loop equations. 
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1 Introduction 



Recently there has renewed interest in the problem of inducing QCD by means of some 
pre-theory. As was proposed by Kazakov and Migdal |1J, such a theory can be potentially 
solvable in the limit of large number of colors, N c , providing the inducing model is that 
of the (self-interacting) matrix scalar field in the adjoint representation of the gauge 
group SU (N c ) on the lattice. The gauge field is attached in the usual way to make the 
model gauge invariant except no kinetic term for the gauge field. The latter circumstance 
differs the Kazakov-Migdal model from the standard lattice Higgs-gauge models. It was 
conjectured in Ref. Jl| that the model undergoes, with decreasing the bare mass of the 
scalar field, a second order phase transition which is associated with continuum QCD 
when the critical point is approached from the strong coupling region. 

To solve the Kazakov-Migdal model in the strong coupling region, Migdal [0] applied 
the Riemann-Hilbert method and derived the master field equation to determine the 
N c — oo solution. An explicit solution of this equation for the quadratic potential is 
found by Gross 0. A surprising property of the master field equation (not yet completely 
understood) is that it admits [§, f§] self-consistent scaling solutions with non-trivial critical 
indices for the non-quadratic potentials. Moreover, the very Riemann-Hilbert method of 
Ref. was developed, in fact, to find such a scaling solution. 

These nice features of the Kazakov-Migdal model are due to the fact that the scalar 
field is in the adjoint representation of the gauge group. For this reason it can be di- 
agonalized by a (local) gauge transformation so that only O(N ) degrees of freedom are 
left and the saddle-point method is applicable as N c — > oo. However, a price for having 
the adjoint-representation field is an extra local symmetry which leads H] at N c = oo 
to local confinement (i.e. the infinite string tension) rather than area law for the Wilson 
loops. 

In the present talk I review the papers || |7|, |S| where the questions of which models 
induce large-iV c QCD with normal area and how to solve these models in the strong 
coupling region were answered^. In Sect. |2] I describe the models, both scalar and fermion 
ones, which induce large-iV c QCD. The mechanism exploited is based on the first order 
phase transition which occurs with decreasing bare mass of the inducing field and is 
associated with freezing the gauge field and the restoration of area law. In Sect. |3] I 
discuss the exact strong coupling solution for the quadratic potential, both in scalar and 
in fermion cases. It is obtained by solving loop equations which turn out to be a useful 
tool for studies of the matrix models. 



2 For a review of other approaches, see recent surveys M nO, and references therein 
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2 The models of induced large- N QCD 



Since the inducing matter field is in the adjoint representation of SU(N C ), Wilson loops 
vanish to each order of the large mass expansion at iV c = oo. This situation is associated 
with local confinement. However, the area law is restored with decreasing bare mass at the 
point of the first order large- iV c phase transition. Its existence can be rigorously proven 
for the single-plaquette adjoint action and is expected for more complicated models on the 
basis of the mean field method. Once the first order phase transition occurs, the proper 
model will induce large- iV c QCD in the continuum. 



2.1 Adjoint scalar model 

A simplest way to induce large- N c QCD is by adjoint representation scalars at Nf flavors 
(i.e. the number of different species). The adjoint scalar model ASM is defined by the 
partition function 



I Watt l MT TT I \ £/=i £. te (-V[*/(*)I+££=i 
'asm = / [[dU^x) [[ [[ d$ f (x)e ' « \ > 

(2.1) 

where the fields $/(x) (/ = 1, . . . ,Nf) take values in the adjoint representation of the 
gauge group SU (N c ) and the link variable U^(x) belongs to the gauge group. The potential 
V[$] is given by 

V[$] = y$ 2 + . . . (2.2) 

where m is the (square of the) bare mass of the scalar field. The original Kazakov-Migdal 
model [J] corresponds to Nf = 1. Notice that the action in Eq. ( |2.1|) is diagonal w.r.t. 
the flavor indices. 

The case of Nf = 1 is a unique one when the matrix $(x) can be reduced to a diagonal 
form at each site of the lattice by means of a gauge transformation. Only at Nf = 1 the 
Riemann-Hilbert method of Ref. is therefore applicable. The alternative method of 
solving ASM with the quadratic potential at strong coupling is based on loop equations 0] 
and can be used at any Nf while gives at Nf = 1 the same result as the Riemann-Hilbert 
method. 



It is worth noting that one can integrate in (|2.1|) over arbitrary hermitean matrices 
rather than over those in the adjoint representation which gives at finite N c a different 
model for the case of a non-quadratic potential. These two modeld should coincide, 
however, as N c — » oo. 
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2.2 Adjoint fermion model 

An alternative to ASM is the adjoint fermion model (AFM) which is defined by the 
partition function 

Nf 

z AFM = /n^nn^/w^/W 6 ^^- ( 2 - 3 ) 

Here Sf[®,^?,U] is the lattice fermion action 

Nf 

S F [V, *,U] = J2 y Z N * tT (Veven(^f(x)^f(x)) 
f=l x 

-j2[^ f (x)p-U,(x)^ f (x + ^Ul(x) +* f (x + riP+Ul(x)* f (x)U> 1 (x)}) (2.4) 

where 

U(M)=mM + ... (2.5) 
is a fermionic analogue of the potential Q2.2Q and m is the bare mass of the fermion field. 
In Eqs. ( |2.3| ), fl2.4|) ^ f{x) is the Grassmann anticommuting N c x N c matrix field while 

I? = r±7 ll (2.6) 

stand for the projectors. The case r = corresponds to chiral fermions while r = 1 is 
associated with Wilson fermions. As is well known, the chiral fermions describe 2 D iVy 
flavors in the naive continuum limit while Wilson fermions are associated with Nf flavors. 



2.3 The induced action 

The above matrix models can be represented in the form of a gauge theory given by the 
partition function 

Z = [ \\dU^x) e' s - Mx)] , (2.7) 

where the induced action for the gauge field U^x) is defined by the integral over <&(x) in 
Eq. fl2~TD or over \&(x) and ^f(x) in Eq. flOj ): 

x f=l 

or 

N f 

e -s ind [£/„(*)] = / JJ JJ ^(^df e ? (2.9) 

X f=i 

with U] given by Eq. ( gg) . 
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For the quadratic potential the result of integrating over or over and ^(x) 
is given by the large mass expansion: 

SMlu]= -^^mi(ir (2 ,o) 



2 y z(r) Vmy 

for scalars or 



, trU(V) I 2 

SW^] = -iV/E '^L/n S P 11^ (2-11) 



/(r)m z ( r ) 



for fermions. In Eq. ( 2.11 ) Sp stands for the trace over the spinor indices of the path- 
ordered product of the projectors (|2.6|) (plus or minus depends on the orientation of the 
link I) along the loop T. 

One easily sees that only single plaquettes survive in the sum over path on the r.h.s.'s 
of Eqs. ( [2,10| ) and ( 2,11| ) if m ~ N 1 /* as N f — > oo, so that the single plaquette adjoint 
action arises in the large- Nf limit: 

S A = -^Y,\^ U ( d P)\ 2 (2-12) 



2 , 



with 



4N f 

(3a = — t for scalars; 

m 4 

(5 A = ^/(l + 2r rj^ for fermions. (2-13) 

This shows of how ASM and AFM induce the single plaquette lattice gauge theory with 
adjoint action. 



2.4 The large-iV phase transition 



The inducing of large- N c QCD relies on the fact [12| that the lattice gauge theory defined 
by the partition function ( |2.7|) with the action ( 2.12[) undergoes the first order large- N c 



phase transition at (3a ~ 2 after which the gauge field U^x) becomes frozen near some 



oo 



mean-field value rj (rj — > 1 as (3 A 

The proof of the existence of the phase transition is based solely on the factorization 
at large-A^c which says that the adjoint action ( 2.12|) is equivalent at N c = oo to the 
single-plaquette fundamental action 



S F [U] = N c pJ2^ tT u ( 9 p) 



(2.14) 



providing the coupling (3 is determined by 



P = (3 A W F (dp-J) 



(2.15) 
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where Wf(Op; 0) stands for the plaquette average 



W F (dp;(3) = 



(2.16) 



Eq. (|2.15 ) can be naively obtained substituting one of two traces in the action fl2.12|) by 
the average. A rigorous proof |L2| is based on the loop equations. 

The existence of the first order phase transition for the action ( |2.12|) with decreasing 
/3a can be seen as follows. Let us solve Eq. ( |2.15| ) for 0(0 a) substituting for (|2.16j ) the 
strong coupling expansion 



5 

W F (dp; 0) = ^ + ^ 



at strong coupling . 



[2.17) 



Eq. ( |2.15| ) possesses at any a trivial solution — 0. However, one more solution appears 
for « 2: 



which matches the weak coupling solution 

~0 - 0A ~ t 



as 0/ 



oo . 



:2.i8i 



;2.i9) 



Notice that -C 1 for the solution ( 2.18|) when 0a ~ 2 so that the strong coupling 
expansion is applicable. 

The adjoint plaquette average 

JU^dU^x) e- s -M^(\trU(dp)\ 2 - 



W A (dp;0 A ) 



JUx^dU^x) e 



which is given due to the factorization by 

W A (dp;0 A ) = (w F (dp-J) 



-,\2 



S A [U] 



(2.20) 



:2.2D 



is depicted in Fig. 1. Since the slope is negative for the solution ( 2.18Q near 0a = 2, 
a first order phase transition must occur with increasing 0a- This negative slop is a 
consequence solely of the positive sign of the second term on the r.h.s. of Eq. ( |2.17| ). The 
predicted value of the critical coupling A , at which the phase transition occurs, obeys 
0\ < 2, as is seen from Fig. 1, to be compared with the result of Monte-Carlo simulations 
0* A = 1.7-1.8. 



2.5 The mean field phase diagram 

At finite Nf the induced actions (|2.10|) or ( p. 11] ) can not be exactly analyzed even at 
N c = oo. An approximate mean field method, which usually works very well in the cases 
of first order phase transitions, was applied to obtain the phase diagram in Refs. H, M . 
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Naively, the mean field approximation consists in substituting the link variable U^x) 
by the mean field value 

[U,{x)f = V ^ (2.22) 

everywhere but one link and writing a self-consistency condition at this link. The self- 
consistency condition is given by the one-link problem 



V = ~ L 2.23 

where 

6a = Jl\ x d$(x)e^ V ^ ^^tr$(0)$(0 + / .) 

/n I #(x)e EiiVtr (^ Wl),+, ' 2E ^ (im ' l) ) 

These naive mean field formulas can be obtained || in the framework of the variational 
method. 

The mean field phase diagram which was obtained by an analysis of Eqs. (|2.23|) and 
( ggg ) is depicted in Fig. 2. At N f = 1 there is no first order phase transition for the 
quadratic potential in the stability region m > 2D. For m < 2D the model is unstable 
and were in the Higgs phase if the stabilizing higher order in $ terms would be added to 
the potential ( |2.2j ). The desired large- N c phase transition appears for Nf > Nf w 30.0 
ASM looks in this region exactly like the single-plaquette adjoint model discussed in the 
previous subsection. 

A similar phase diagram for AFM is depicted in Fig. 3. Now there is no Higgs phase (or 
an unstability region for the quadratic potential) due to the fermionic nature of inducing 
fields. For the cases of chiral and Kogut-Susskind fermions the first order phase transition 
is present already for Nf = 1 while the result for Wilson fermions is less certain. 



2.6 Area law versus local confinement 

At the point of the first order large- N c phase transition, the area law behavior of the 
(adjoint) Wilson loops which is associated with normal confinement is restored in ASM 
or AFM analogously to the single-plaquette adjoint action JL^ . 



In order to see this, let us consider the adjoint Wilson loop which is defined by 

W A (C) = /^(\trU(C)\ 2 -l)\ (2.25) 



where the average is understood w.r.t. the same measure as in Eq. ( |2.1| ) or in Eq. ( |2.3| ). 
Alternatively, one can average w.r.t. the induced actions ( 2.1(J| ) ( [2.11| ) which recovers at 



J Such a phase diagram is compatible with the Monte-Carlo studies of Rcf. [fL3| 
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Nf = oo the single plaquette adjoint action (|2.12 ). In this limiting case the following 
extension of the factorization formula ( |2.21 ) holds at N c 



oo: 



W A {C-(3 A ) = (W F (C-J) 



(2.26) 



where Wf(C; (3) is defined by the same formula as ( 2.16|) with dp replaced by an arbitrary 
contour C and f3 versus [3a given by Eq. ( |2.15p . 

Since (3 = for (3a < (3% W A (C) vanishes in this region due to Eq. ( |2.26| ) except the 
loops with vanishing minimal area A m i n (C): 



[2.27) 




W A (C) = 5 0Amin( c) + O 

On the contrary, the area law with the string tension 

K A ((3 A ) =2K F (P(J3 A )) 



(2.28) 



holds for f3 A > (3* A when Eq. (|2.15|) possesses the non-trivial solution. An extension of 
these formulas to finite Nf is given in Ref. 0. 

While the first order phase transition associated with the restoration of area law looks 
similar for ASM and AFM, the continuum limits should be approached in different ways. 
For ASM the continuum QCD is reached at the line of second order phase transitions 
which separates the area law and Higgs phases provided that one approaches it from the 
area law phase. For AFM there is no Higgs phase and continuum QCD is reached as 
0. 



in 



3 Loop equations at strong coupling 

The loop equations of ASM or AFM relate the closed adjoint Wilson loop ( |2.25| ) to 
the open ones with scalars or fermions at the ends. The loop equations are drastically 
simplified at N c = oo in the strong coupling region where the closed loops obey Eq. ( [2. 27 ). 
The exact solution can be obtained in both cases for the quadratic potential when the loop 
equations turns out to be equivalent to those for the hermitean and complex one-matrix 
models, respectively. 

3.1 Loop equations for arbitrary potential 

The generic object which appear in the loop equations are open Wilson loops 

5 ff G x (C xy ) = t±r tr Uf{x)U{C xy ) x _l U\C xy ) \ \ . (3.1) 
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The appearance of the <5-symbol w.r.t. the flavor indices / and /' is due to the fact that 
the interaction terms in the action entering Eq. ( |2.1| ) are diagonal over the flavor indices. 

The loop equations of ASM result from the invariance of the measure in Eq. ( |2.1| ) 
under an arbitrary shift of and reads 



D I 1 1 

- £ ( y c tr ( $/ ^ + v)u(c {x+IJ , )x c xy ) — ^-^ u\c {x+ll)x c xy 



xy. 



H=—D \-"C 



(itr (^^—L^tr (—L^t^)^ (3.2) 

where the path C^^jjC^j on the l.h.s. is obtained by attaching the link (x, fi) to the path 
C xx at the end point x as is depicted in Fig. 4. I have omitted additional contact terms 
which arise at finite N c due to the fact that $ belongs to the adjoint representation, so 
that Eq. ( |3.2|) is written for the hermitean matrices. This difference disappears, however, 
as N c — > oo. 

The analogues of Eqs. Q and (p]) for AFM read 

SffG'HCJ = (lt,(^( C J A ,_^ fM g t(C,))) (3.3) 
where z and j are spinor indices, and 

tr (VfW^mxWMmC.,) x _ *^ f> {y) U\C xy ) 

-lt(w tT ( P t*f( x + nme^Csy) _ U\C {x+ , )x C xy 



+p^ f (x - #(c Nl g A _ f f ll)l f{y) u\c^ )x c xy ))) 



■^c v v y A 2 - ^ f {y)^ f{y)> N c ^ - * f (y)* f (y) 
The matrix multiplication over the spinor indices is implied in this equation. 



3.2 Loop equations at large N c 



The path C xy on the r.h.s. of Eq. fl3.2|) (or Eq. ( [3^4] )) is always closed due to the presence 
of the delta-function. The explicit equation for the case of vanishing (or contractable) 
contour C xx = at large N c , when the factorization holds, reads 

du V'(uj) ^ n ^ ^ r2 



Ci 2-ni (A — oj) 



E w - 2DG X (1) = E{ (3.5) 



9 



where 

M^ tr (^7M)K (Gj(0) + 1) (3 ' 6) 

with G\ is defined by Eq. (|3.1| ). I have denoted the one-link average by 

G A (1) = G X {C (X ^ )X ) (3.7) 

since the r.h.s. does not depend on x and \x due to the invariance under translations by 
a multiple of the lattice spacing and/or rotations by a multiple of 7r/2 on the lattice. 
The contour C\ encircles singularities of E w so that the integration over u on the l.h.s. of 
Eq. ( |3.5|) plays the role of a projector picking up negative powers of A. 

For C xx 7^ 0, the averages of a new kind arise on the r.h.s. of Eq. (|3.2|) (or Eq. ( |3.4fj ). 
However, these averages obey at N c = oo the following analogue of Eq. fl2.27| ) 

(3.8) 

i.e. vanish for C xx ^ 0. 

Hence, the loop equation for C xy ^ at N c = oo reads 

/_L tr (y'(<f> f (x))U(C xy ) — I {y) rf( c J)) ~ E G x (C ix+ , )x C xy ) = (3.9) 

independently of whether C xy is closed or open. Therefore, the r.h.s. of the loop equation 
in nonvanishing at N c = oo only for C xy = (modulo backtrackings) when the proper 
equation is given by Eq. (|3.5|) . 

Finally, the fermionic analogues of Eqs. (|3.8|) and ( p.6|) read 

A \\ 



N c tT { U{Cxx) \*-* f (x)y f (x)) N c tT { Ui ( Cxx) \2-* f (x)* f (x 

\N 2 



S 0Amm(c) El + o(^) (3.10) 



and 

^(^( y-W) ))'^" 1 " (3 ' u) 

with G\ defined by Eq. ( p.3|) . Therefore, the r.h.s. of Eq. (|3.4j ) involves at N c = oo only 
E\ similarly to the scalar case. 



3.3 The quadratic potential 

The quadratic potential is always solvable, even in the non-diagonizable cases, for the 
following reasons. Let us consider the one-link correlator 

,1 JdU e N '^i*>" x ' u, ±-tTrU X rtf 
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where t a (a = 1, . . . , N% — 1) stand for generators of SU(N C ) which are normalized by 

1 



i- tr t a t b = 5 ab . (3.13) 

Nr. ^ ' 



At N r = oo the formula 



i- tr t a Uxf>U^) u = tr f ^ , (3.14) 

where A is a constant to be determined below, can be proven for (p and \ given by the 
master field for the quadratic potential analyzing the large mass expansion. The point is 
that terms like tr t a (p k j, with k > 1 never appear for the quadratic potential. Analogously 
it can be proven that 

-^-ti t a U— — U j ) =A-Ltrt a — — . (3.15) 

For Nf = 1 Eqs. ( |3.14 ) and ( 3.1 5| ) recovers the ones of Ref. [14 . 
For G x (C xy ) defined by Eq. Q), Eq. flQ5|) implies 

G{C xy ) = A L E\ (3.16) 

where L is the algebraic length (i.e. the one after contracting the backtrackings) of C xy . 
The fermionic analogue of this formula reads 

Gt{C xy ) = A L E X ( I] P±) (3.17) 

VeC„ s / 

where the plus or minus signs correspond to the direction of the link I which belongs to 
the contour C xy . The spin factor will provide below the cancellation of the projectors in 
Eq. (0). 

The constant A can be determined by substituting the ansatz ( |3.16|) into the C xy ^ 
loop equation (|3.9| ) which simplifies for the quadratic potential as 

D 

mG\(C X y) — ^2 G\(C( x+ n) x C xy ) = 0. (3.18) 

fi=-n 

The ansatz (|3.16|) satisfies this equation for any C xy ^ providing 

A = 2 = (3.19) 

m + Jm 2 + 4(1 - 2D)a 



where a = 1 for scalars or 



a = P+P- = r 2 - 1 (3.20) 
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for fermions (a = — 1 for chiral fermions and a = for Wilson fermions). 

The remaining function E\ can now be determined from Eq. ( |3.5|) which for the 
quadratic potential reads 

rhE x = E 2 X , rh = m- 2DA (3.21) 

and coincides with the loop equation for the Gaussian hermitean one-matrix model (for 
a review, see Ref. |HJ and references therein). The solution of Eq. (|3.21 ) which satisfies 



E x -> \ as A -> oo , (3.22) 
A 



as it should be due to the definition (|3.6|), is unambiguous 



2E X = m\- rh^/X 2 - 4/m . (3.23) 

The imaginary part 



SEx = p(A) = —mJA/rh - A 2 for - 2/Vm < A < 2/Vm (3.24) 

4-7T V 

recovers the solution [[|. Analogously, the C xx = loop equation for AFM with the 
quadratic potential is reduced to the loop equation for the complex one- matrix model |15|] . 

One should not be surprised that the exact strong coupling solution for the quadratic 
potential does not depend on Nf which is a consequence of the peculiar behavior of the 
Wilson loops ( |2.27| ). This independence does not contradict to the fact that the first 
order phase transition discussed in Sect. |^ occurs only for Nf > Nj. The point is that 
the strong coupling solution is not sensitive to the phase transition which occurs due to 
another (thermodynamic) reason. 



3.4 Interpretation as the ID tree problem 

A question arises what combinatorial problem are the exact solutions of the previous 
subsection associated with? To answer, let us consider the open loop correlator 

6ff>G(C xy ) = (i-tr (® f {x)U{C xy )$ fl {y)U\C xy j) ) . (3.25) 

which is nothing but the A -2 term of the expansion of (3.1) in A~ x . At = oo the 
standard sum-over-path representation of G{C xy ) reads 

/ 2 \'( r )+i 

G(C xy ) = E - W A (C xy T yx ) (3.26) 

1 yx 

where the contour T yx forms together with C cy a closed loop passing x and y. 
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Since in our case the formula ( |2.27| ) associated with the infinite string tension holds, 
Y yx must coincide with (C^) -1 (i.e. passed in opposite direction) modulo backtrackings 
of T. Eq. ( |3T26D then yields 

/ 2 \'( r ) +1 

G«W = E(-) (3^) 

1 yx 

where the sum goes over contours of the type depicted in Fig. 5.0 
The fermonic analogues of Eqs. ( |3.25| ) to ( |3.27|) read 



i.i 



i J(r)+i 

G ij (C xy ) =£(-) W A (C xy T yx ) ( J] Pt ] (3-28) 



and 



m 



J(T)+1 



G%)=eQ) Spinal (3-29) 

with summing again over contours depicted in Fig. 5. 

The proper combinatorial problem is, therefore, that of summing over 1-dimensional 
trees embedded in a D-dimensional space. The exact solution of loop equations for the 
quadratic potential represents the solution to this problem: 

2D — 1 

G(C xy ) = A L , (3.30) 

m(D - 1) + D^Jm 2 + 4(1 - 2D)a 

where L is the algebraic length of C xy and A is defined by Eq. (|3.19| ) . Such a dependence 
on A is evident from the representation ( 3-27 ) (or ( |3.29|) ) since the trees are uniformly 
distributed along C xy . 

As is already mentioned in the previous section, this solution coincides for scalars 
with that of Ref. ||. Eq. ( |3.30| ) takes an especially simple form for Wilson fermions 
when the backtracking parameter a, given by Eq. ( |3.20|) , vanishes so that there are no 
backtrackings. For chiral fermions when a = — 1 the solution ( |3.30j ) coincides with that 
of Ref. [|T6[ for the case of the fundamental representation and vanishing constant in front 
of the plaquette term in the action. The point is that Wilson loops vanish in this case as 
well (except for those with vanishing minimal area) and exactly the same combinatorial 
problem of summing the diagrams of the type depicted in Fig. 5 emerges. 
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Figures 



W A (dp;/3 A ) 




Fig. 1 The two solutions of Eq. (2.15). The line which starts at [3a = 2 is associated with the 
solution ( 2.1 8| ) . Since the slope is negative for this solution near (5a = 2 (and Wa *C 1), 
the first order phase transition occurs at some (3* A < 2 so that the actual behavior of 
Wa(9p;/3a) is depicted by the bold line. 
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Fig. 2 The mean field prediction for the phase diagram of ASM. The bold lines which 
bounds the local confinement phase is that of first order phase transitions. The line which 
separates the area law and Higgs phases is that of second order phase transitions. 
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3 The mean field prediction for the phase diagram of AFM. The bold line which separates 
the local confinement and the area law phases is that of first order phase transitions. 
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Fig. 4 The graphic representation for G\{C xy ) (a) and G\(C( x+ ^ x C xy ) (b) entering Eq. (J37 
The empty circles represent ^/(x) or <frf(x + /x) while the filled ones represent x-^ ,(y) • 
The oriented solid lines represent the path-ordered products U(C xy ) and U(Cr a 
The color indices are contracted according to the arrows. 
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Fig. 5 The typical paths T yx which contribute the sum on the r.h.s. of Eq. fljDTp (and 
Eq. ( |3.29D ). These coincide with C xy passed backward modulo backtrackings which 
form a ID tree. The sum over T yx is reduced to summing over the backtrackings. 
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